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We describe the theoretical advances that influenced the experimental creation of 
vibrationally and translationally cold polar 40 K 87 Rb molecules [TJ [2] . Cold molecules 
were created from very-weakly bound molecules formed by magnetic field sweeps 
near a Feshbach resonance in collisions of ultra-cold 40 K and 87 Rb atoms. Our 
analysis include the multi-channel bound-state calculations of the hyperfine and 
Zeeman mixed X 4 £ + and a 3 S + vibrational levels. We find excellent agreement with 
the hyperfine structure observed in experimental data. In addition, we studied the 
spin-orbit mixing in the intermediate state of the Raman transition. This allowed 
us to investigate its effect on the vibrationally-averaged transition dipole moment to 
the lowest ro- vibrational level of the X 1 S + state. Finally, we obtained an estimate of 
the polarizability of the initial and final ro- vibrational states of the Raman transition 
near frequencies relevant for optical trapping of the molecules. 



I. INTRODUCTION 

The recent successful creation of a high phase-space-density gas of polar 40 K 87 Rb 
molecules [TJ |2] has been based on both new experimental and theoretical advances in 
manipulating and understanding properties of such molecules. This opens up the possibility 
of studying collective phenomena that rely on the long-range interactions between polar 
molecules. Future experiments can be envisioned in both weakly confining optical traps as 
well as optical lattices. 

Our goal in this paper is to describe some of the theoretical advances that influenced the 
experimental creation of vibrationally and translationally cold polar 40 K 87 Rb molecules. In 
particular, we theoretically analyze various factors that can affect this creation including 
the multi-channel description of the initial, intermediate, and final states of the formation 
by Raman transitions. 

In a previous paper [3] we made the first steps towards obtaining practical guidelines for 
photoassociatively producing low v vibrational states of heteronuclear KRb. We calculated 
the electronic transition dipole moments between ground state of the KRb molecule and 
excited states. In addition, we obtained the permanent dipole moments of the polar X 4 S + 
and a 3 E + ground states. A relativistic electronic structure code was used. 

In a second paper [I] we discussed the possibility of creating X 4 S + molecules starting 
from doubly spin-polarized K and Rb atoms via two-photon photoassociation. We assumed 
that colliding atoms are initially in the doubly spin-polarized state, which only allows them 
to bond in the ground configuration a 3 S + potential. We then found that there is a viable 
route from the doubly spin-polarized colliding atoms to the vibrationally ground X 4 S state 
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FIG. 1: The ground and excited state potential energy curves of KRb that are used for multi- 
channel modeling as a function of internuclear separation R. Here ao is the Bohr radius of 0.0529 
nm. 

via the excited A 1 E + and b 3 II states, when they are mixed through spin-orbit interactions. 
This mixing allows us to describe these potential by the Hund's case (c) coupling scheme as 
2(0 + ) and 3(0 + ) states and calculate the Raman transition rate from the triplet state to the 
excited Q = + state followed by a downward transition to the ground singlet state. This 
process is absent in homonuclear dimers, since the additional gerade-ungerade symmetry 
prevents it. In the notation m(fl) the number in parenthesis is the projection Q of the total 
electronic angular momentum on the internuclear axis and the number in front labels the 
order of states with the same Q. The role of black-body radiation in redistributing population 
among ro- vibrational levels of the singlet X X E + and triplet a 3 £ + states was investigated. 

Here we further search for an efficient production mechanism using a multi-channel de- 
scription of both ground and excited states. We assume that KRb molecules are initially 
in the weakly-bound near-threshold vibrational states formed by a magnetic Feshbach res- 
onance in collisions between ultracold 40 K and 87 Rb atoms. In our coupled-channel calcu- 
lation of the ground state ro-vibrational structure we used the most accurate ground state 
potentials available from Ref. [5] . In section [Tf] we analyze this structure, and perform a 
comparison with available high-precision measurements [U [2] . 

Vibrationally cold molecules are preferably made by transfering population from a Fes- 
hbach molecular state. For this transfer we selected the pathway that has been proposed 
by Stwalley [6], which forms vibrationally cold KRb molecules starting from the highly ex- 
cited vibrational states using one optical Raman transition and intermediate vibrational 
levels of the 3(1) potential. This mechanism was also used to create vibrationally cold 
RbCs molecules in Ref. [7]. Reference [8] has reported an analysis of perturbations of the 
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vibrational levels of the 3(1) potential due to spin-orbit interactions with the neighboring 
potentials in RbCs. 

For the vibrational levels used as intermediate states the 3(1) potential can to first order 
be described as the nonrelativistic 2 3 S + state. More accurately the strong non-adiabatic 
interaction with the neighboring 2(1) and 4(1) potentials has to be taken into account. 
Alternately, we can view this coupling as being due to the spin-orbi t in teraction between 
the nonrelativistic 2 3 £ + , l 3 n, and l 1 !! potentials. Therefore in Section III we perform multi- 
channel calculations of the ro- vibrational structure and the vibrationally-averaged transition 
dipole moments to the ground state levels. The three intermediate excited states that are 
of interest have Q — 1 symmetry and are shown in Fig. [T] together with the ground state 
potentials of KRb. In addition, panel a of Fig. [I] shows the pathway to form v — a 3 S + 
molecules starting from a gas of Feshbach molecules. Similarly, panel b shows the pathway 
to v — X X S + molecules. 

The excited potentials in Fig. [l] were constructed from RKR data [9j [TO] and as well as 
from our ab initio calculations. The two attractive 3(1) and 4(1) potentials dissociate to the 
K(4s)+Rb(5p3/2) atomic limit whereas the 2(1) potential dissociates to the K(4s)+Rb(5pi/2) 
limit. At short internuclear separations the potentials can be approximately described by 
the Hund's case (a) n +1 A symmetry, where A is the projection of the electron orbital 
angular momentum along the internuclear axis and S is the total electron spin. At longer R 
relativistic effects are important, where the curves can only be described with Hund's case 
(c) m(Q = 1) labeling. 

In Section |III| our multi-channel calculation is, however, based on the nonrelativistic 
excited potentials 1 3 II, 1 II, and 2 3 S, electronic transition dipole moments to the ground 
states, and the spin-orbit coupling matrix elements. Some RKR data for the more deeply 
bound vibrational levels of the l 1 !! and 2 3 E are available p.110]. We extend this information 
by ab initio electronic structure data and by known long range dispersion coefficients from 
Ref. [IT]. 

The relativistic configuration interaction molecular orbital restricted active space (MOL- 
RAS-CI) method has been used to calculate potential energy curves, permanent and transi- 
tion electric dipole moments of the KRb heteronuclear molecule as a function of internuclear 
separation. We combine this calculation with multi-channel ro-vibrational structure calcu- 
lation to obtain Frank-Condon factors between the ground and excited states of KRb. 

Finally, in Section IV we describe our calculation of the polarizability of vibrational levels 
of the X 1 S + and a 3 S + states. 



II. COUPLED-CHANNEL CALCULATION OF THE GROUND STATES 

In the Raman transition the initial and final bound vibrational levels belong to the ground 
X : S + and a 3 E + states. In KRb these states dissociate to the same [Ar]4s( 2 S) + [Kr]5s( 2 S) 
atomic limit. The two states are coupled via hyperfine interactions: the Fermi-contact and 
electron and nuclear Zeeman interactions for each of the constituent atoms. For each atom 
the Fermi-contact interaction couples its electron spin, here 1/2, to its nuclear spin. The 
Zeeman interaction is non-zero since an external magnetic field is used to create Feshbach 
KRb molecules. The corresponding Hamiltonian for vibrational states in this coupled system 
has been discussed in Refs. [121 [13]. For this paper we include the effect of the weak magnetic 
dipole-dipole interaction perturbatively but do not include the effect of the second-order 
spin-orbit interaction |14j . 
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Numerically solving for eigen pairs of this system is called a coupled-channel calculation. 
For our calculations we used the electronic potentials of Ref . [5] . The atomic masses for 40 K 
and 87 Rb are taken from Ref. [TS], the Fermi-contact term values and electronic g-factors are 
from Ref. [TB] and, finally, the nuclear magnetic moments are from Ref. [TTJ . The nuclear 
spin of 40 K is 4 and that of 87 Rb is 3/2. 

The vibrational wave functions can be labeled by three nearly conserved angular mo- 
mentum quantum numbers. These are the relative orbital angular momentum £ of the two 
atoms, its projection mt> along the magnetic field direction, and the projection Mp of the 
summed atomic angular momentum F = f a + f b along the same direction. Here, f a and /& 
are total angular momenta of each atom. In fact, we can also write F = 8 + 1, where S and 
/ are the total electron and nuclear spin, respectively. 

For our calculations we neglect coupling between states with different £miMp. In the 
absence of the magnetic dipole interaction the 2£ + \ levels with the same £ and Mp quantum 
numbers are degenerate. For ultracold atoms and molecules we can limit ourselves to £ = 
and 2. As the Feshbach molecule created in [21 [18] has £ = and Mp = — 7/2 we will limit 
ourselves to Mp = —11/2 to —3/2 as in a Raman transition Mp can change by upto two 
units depending on the polarization of the light beams. 

Figure [2] shows the eigen energies and properties of eigen functions of a coupled-channel 
calculation with Mp = —7/2 and £ = at a magnetic field of 5=545.9 G near the dissocia- 
tion limit. The figure makes evident that mixing between the singlet X X E + and triplet a 3 S + 
states becomes strong within 150 GHz of the limit. The vertical axis shows the fraction of 
the a 3 E + character in the wave function. A value close to zero (one) corresponds to a state 
primarily described as a X 1 S + (a 3 £ + ) level. Our calculations finds non-degenerate hyper- 
fine and Zeeman structure, which are grouped in Fig. [2] by a vibrational quantum number. 
Within each group there are twelve sublevels, which for weak mixing corresponds to three 
singlet and nine triplet sublevels. The vibrational quantum number is labeled according to 
the vibrational quantum number of the uncoupled singlet and triplet potential. The v——l 
level corresponds to the last uncoupled bound level. The grouping is valid in this case 
because the two potentials have the same long-range dispersion potential and sufficiently 
similar scattering lengths that the spacing between vibrational levels of the two potentials 
is nearly the same in the range of energy shown. For binding energies less than 3 GHz 
the levels with different vibrational quantum numbers intermix. The level relevant for the 
experiment of [2] , which is the initial state of both Raman transitions in Fig. [TJ is the most 
weakly bound sublevel with Mp = —7/2 and has a 0.23 MHz binding energy at 5=545.9 G. 
Here, this state has 80 % a 3 S + character. The weakly bound levels are now fully discussed 
in Ref. [IS]. 

The final states of the Raman transitions are v — levels of the X 1 S + and a 3 S + potentials. 
Figure [3] shows the rotational hyperfine and Zeeman structure of the v — level of the a 3 S + 
potential at 5=545.9 G as calculated with the coupled-channel method. The black lines are 
the £=0 bound states and the red lines are the £=2 bound state. Notice though that the 
rotational energy splitting is smaller than that due to the hyperfine and Zeeman interaction. 
This leads to overlapping spectral features. The hyperfine structure of the two partial waves 
is nearly identical. The main difference is that each £ = 2 hyperfine feature has within it 
three lines, which are not resolved in the figure. The splitting between these lines is less 
then 0.3 GHz and is due to the magnetic spin-spin dipole interaction, which partially lifts 
the me degeneracy of the projection quantum number of £. 

The experiments of Refs. [UI2] have located more than ten sublevels of the v=0 vibrational 
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FIG. 2: Fraction of the a 3 S + character of the coupled-channel eigen states of the 40 K 87 Rb molecule 
as a function of their energy at a magnetic field of 5=545.9 G. Bound states with Mp = —7/2 
and t = are shown. For energy between — 300GHz and — 15GHz the levels can be grouped 
by vibrational quantum number v of the uncoupled X 1 S + and a 3 S + potentials with v = —1 
corresponding to the most weakly bound levels. Panel a displays vibrational levels from v = —4 
to v = —9 and panel b is blowup of the near threshold region. Zero energy corresponds to the 
dissociation energy with both 40 K and 87 Rb in the energetically lowest hyperfine state. 

level of the a 3 E + potential. We have compared the calculated hyperfine structure of the v =0 
level of the triplet state, shifted up by +15.0 GHz, with the experimental energies of Ref. [2J, 
marked by the crosses and triangles in Fig. [3| The agreement is good. 

Figure [4] shows the hyperfine and Zeeman structure of the v=0 vibrational state of the 
ground X X S+ state of 40 K 87 Rb at 5=545.9 G shifted up to +0.4014 GHz such that the 
energetically lowest i = Mp = —7/2 level coincides with the experimental data indicated 
by triangles in Fig. |4| The two panel show energy levels for the lowest two even partial 
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FIG. 3: The hyperfine and Zeeman structure at 5=545.9 G of the £=0 (black lines) and 2 (red lines) 
rotational levels of the v = vibrational state of the a 3 S + potential of 40 K 87 Rb. The crosses and 
triangles indicate the experimentally observed energies from Ref. [2] . The theoretical energies have 
been shifted up by +15.0 GHz to coincide with the experimental data. Zero energy corresponds 
to the dissociation energy of both K and Rb in the energetically lowest hyperfine state. The 
levels are grouped by the projection quantum number Mp. Each £ = 2 hyperfine feature contains 
three lines, which on the scale of the figure are barely resolved. The splitting is on the order of 0.1 
GHz and due to the magnetic spin-spin dipole interaction, which partially lifts the rri£ degeneracy 
of the projection quantum number of £. 



waves. The singlet potential has to first order no hyperfine structure due to the Fermi 
contact interaction and electronic Zeeman interaction. Hence the structure in Fig. [4] is 
predominantly due to the nuclear Zeeman interaction of both atoms and is on the order of 
a few MHz at 5=545.9 G. In fact, the nuclear Zeeman energy for the X X S + state is given 
by 

Ez = -{gi,Km K + gi,Rbm Rb )fi N B , (1) 

where mx and are the projections of the nuclear spin of 40 K and 87 Rb along the magnetic 
field direction, respectively. We then have M F = rriK + mm- The gj t K and gi^m are nuclear 
g-factors and fijy is the nuclear magneton. The g-factor of K is about a factor of two smaller 
than that of Rb. In Fig. [4] each line can be labeled by and m^. The smaller splittings 
between lines correspond to levels with different rnx for the same m^. The larger gaps 
between groups of levels correspond to different miu>- Each line of Fig. [3], panel b, contains 
five unresolved components corresponding to the five 1=2 sublevels. Unlike the 1=2 lines of 
the a 3 E + state, the spin-spin dipole interaction here is zero to the first order and the second 
order contribution is very small. 

For the X 1 E + potential Ref. [2] has observed a single hyperfine component for the £=0 and 
2 rotational state of the v =0 vibrational level. The experiment and theory agree to ~ 401.4 
MHz in an energy levels. The theoretical energy difference between £=0 and £=2 components 
of the M F = —7/2 levels is 6.68376 GHz, which agrees well with the experimental value of 
6.6836(5) GHz given in [2]. 
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FIG. 4: The hyperfine and Zeeman structure of the v =0 £=0 level (panel a) and 1=1 (panel b) of 
the X 1 ^ state of 40 K 87 Rb at B = 545.9 G. The triangles with one-standard-deviation error bar 
indicate the experimentally observed energies from Ref. [2j. The theoretical energies have been 
shifted up by +0.4014 GHz such that the energetically lowest t = Mp = —7/2 level coincides 
with the experimental data. The energy regions shown in the two panels are not the same. Zero 
energy corresponds to the dissociation energy of both 40 K and 87 Rb in the energetically lowest 
hyperfine state. The levels are grouped by the projection quantum number Mp. 

III. MULTI-CHANNEL CALCULATION OF THE EXCITED STATES 

In this section we model the ro-vibrational motion of the 40 K 87 Rb molecule in excited 
electronic potentials, which are used as intermediates to create vibrationally cold molecules 
[2] . In particular, we focus on the need to use a multi-channel description of the vibrational 
structure that includes coupling between the electronic potentials. The origin of coupling can 
be explained from either the relativistic spin-orbit interaction, which couples non-relativistic 
2S + 1 A ± Born-Oppenheimer (BO) potentials [20] or non-adiabatic mixing of relativistic $7 ± 
potentials. Some of the relativistic potentials as function of internuclear separation R are 
shown in Fig. [T] However, in this Section we first calculate multi-channel vibrational energies 
based on the non-relativistic potentials. We then discuss the effects of the multi-channel 
calculation on the vibrationally-averaged transition dipole moments. Due to the already 
complex nature of these calculations, we have not included the contributions of hyperfine, 
Zeeman, or coriolis interactions. 

Note that as stated before for a non-relativistic potential the quantum number S, cor- 
responding to the total electron spin S, and A, corresponding the absolute value of the 
projection of the total electron orbital angular momentum, is conserved. For a relativistic 
calculation only Q, the absolute value of the projection of the summed electronic orbital 
and spin angular momentum, are conserved. The ± superscript, only relevant for A = 
or Q = states, distinguishing states with opposite reflection symmetries. For a given Q 
symmetry nonrelativistic potentials that satisfy — S <Q — A < S, are coupled by spin-orbit 
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FIG. 5: Panel a: the -R-dependent off-diagonal spin-orbit matrix elements between the 2 3 £ + 
and l 1 !! potentials (solid line), the 2 3 £ + and 1 3 II potentials (dashed line) and l 1 ]! and 1 3 II 
potentials (dash-dotted line). For the spin-orbit interaction between f2 = 1 potentials diagonal 
matrix elements are zero. Panel b: the independent transition dipole moment between the X 1 S + 
and two 1 II states. The dipole moments are in units of eao, where e is the electron charge. 



interactions. Strong mixing occurs when the energy splitting between 25 + 1 A ± potentials is 
on the order of the spin-orbit interaction energy. 

In our mult i- channel calculation we use the X 1 S + non-relativistic ground-state BO po- 
tential from Ref. [5] and the 2 3 £ + , l 1 !!, and 1 3 II excited potentials from Ref. [2Tj . All three 
excited potentials dissociate to the K( 2 S)+Rb( 2 P) limit. The previously unknown spin-orbit 
coupling matrix elements and electronic dipole moments are obtained from our MOL-RAS- 
CI calculations. The panel a in Fig. [5] shows the independent spin-orbit matrix elements 
between the 2 3 £ + , 2 1 IT, and 1 3 II states. For the spin-orbit interaction between Q=l po- 
tentials diagonal matrix elements are zero. As we will be interested in vibrational levels 
near the bottom of the 3(1) potential spin-orbit coupling to 2 1 IT potential can be neglected. 
At large R the matrix elements approach A/3, where A is the spin-orbit splitting of the 
2 P state of 87 Rb. The panel b in Fig. [5] shows the independent transition dipole moment 
between the X 1 S + and the l 1 !! and 2 1 IT states. The dipole moments between the singlet 
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FIG. 6: The fraction of 3 n (blue), 3 £ + (red), and 1 I1 (green) character of the multi-channel = 1 
eigen states of the 40 K 87 Rb molecule as a function of energy. The level indicated by the label v = 23 
corresponds to the v=23 bound state of the 3(1) potential used in [2]. Zero energy corresponds to 
the dissociation energy of both 40 K and 87 Rb in the energetically lowest hyperfine state. Hence 
the energy can be interpreted as the photon energy needed to make upward transition shown in 
the right panel of Fig. [T] 

X : S + and the triplet 2 3 S + and 1 3 II states are strictly zero. 

Figure [6] shows the results of a multi-channel calculation of the vibrational levels of 
the Q — 1 excited states. We have coupled three channels, 2 3 £ + , 2 1 I1, and 1 3 IT, and, 
in addition, included the rotational potential h 2 J(J + l)/{2jj,R 2 ) with rotational quantum 
number J = Q — 1 to each channel. The vibrational energy on the horizontal axis is relative 
to the dissociation energy of ground state 40 K and 87 Rb and thus corresponds to the photon 
energy for the upward part of the Raman transition. Each vibrational level is represented 
by three circles, of different color but positioned at the same eigen energy, corresponding 
to the fraction of the three-channel wavefunction that is in the 2 3 S + , 2 1 I1, and 1 3 I1 state, 
respectively. The sum of these fractions adds up to one. A level with a fraction one in a 
single channel corresponds to a vibrational level of an unperturbed 2S+1 A ± potential. 

The energy range in Fig. [6] spans from just below the bottom of the 3(1) potential (See 
Fig. I]) to the atomic K( 2 S)+Rb( 2 P3/2) limit. The v = vibrational level of the 3(1) state 
can Be identified at 9200 cm -1 . For energies E > 9200 cm -1 levels with a large fraction of 
3 S + character appear. For E < 9200 cm -1 only vibrational levels of the 2(1) or 3 I1 state 
exist. Similarly at E w 10900 cm -1 eigen states with a large fraction in the 1 I1 state appear. 
This corresponds to the bottom of the 4(1) potential. Interestingly, for all eigen states with 
energy E < 10900 cm -1 the levels have a small amount, < 0.01, of 1 II character. The 2(1) 
and 3(1) states have a small 1 I1 admixture due to second-order spin-orbit mixing. 

In other energy regions the characterization of levels is less clear. For example levels 
with energy larger than E ~ 10200 cm^ 1 have non-negligible contributions from the 2 3 S + 
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FIG. 7: Transition dipole moment from the v=0, J=0 X 1 !]" 1 " level to the excited 0=1 levels of 
the 40 K 87 Rb molecule as a function of the excited-state energy. Zero energy corresponds to the 
dissociation energy of both 40 K and 87 Rb in the energetically lowest hyperfine state. The bound 
v=23 level of the 3(1) potential, used as intermediate state in [2], is marked. The calculated dipole 
moment to the v=23 level is 0.018 a.u. 

and 1 3 II states. These BO potentials of Ref. (21] cross at this energy and the spin-orbit 
interaction mixes the two symmetries. The intermediate vibrational level used in Ref. (2] 
and indicated by v = 23 in Fig. [6] is such a mixed state. From our calculation we find that 
it has a 79% 3 S+, 21% 3 IT, and 0.2% x n character. The closeness of the v = 23 level to the 
avoided crossing and the theoretical uncertainties in its location make the precise fractions 
uncertain. In fact, a few hundred cm -1 upward shift could potentially remove all 3 S + and 
3 I1 mixing and the v = 23 level becomes a nearly pure 2 3 S + vibrational level. The 1 I1 
character, however, is not expected to change significantly 

Figure [7] shows the transition dipole moments between vibrational levels of the multi- 
channel Q — 1 calculation described above and the v = 0J = £ = ro- vibrational level of the 
X : S + potential assuming the electronic dipole moment shown in Fig. [H] This vibrationally- 
averaged dipole moment describes the downward part of the Raman transition. The energy 
region is as in Fig. [6] and the changing character of the excited state vibrational levels is 
reflected in the transition dipole moments. The vibrational averaged dipole moment from 
the singlet X X S + state is only nonzero if the multi-channel vibrational levels contains 1 I1 
character. A larger character leads to a larger transition dipole moment. The start of 
the vibrational series of the 3(1) and 4(1) potential are clearly visible in the figure. The 
intermediate vibrational level used in Ref. [2] is again indicated by v = 23. From our 
calculation we find a dipole moment of 0.018 a.u. for the transition from this level. 
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IV. DYNAMIC POL ARIZ ABILITY OF THE GROUND STATE VIBRATIONAL 

LEVELS 

We examine the dynamic polarizability a of the KRb molecule as a function of laser 
frequency, lu, and ro-vibrational quantum numbers, v, J, of the uncoupled ground state 
a 3 E + and X 1 E + potentials. The real part of the dynamic polarizability, among other things, 
determines the depth of the trapping potential seen by a molecule as 

V = —Re(a(huj,v)) x J, (2) 

where I is the intensity of the laser fields at frequency u. The imaginary part of a describes 
the spontaneous or any other decay mechanism that leads to loss of molecules from the trap. 
Based on its knowledge, laser frequencies can be selected to minimize decoherence effects 
from loss of molecules due to spontaneous or laser-induced transitions. 

The dynamic polarizibility of a ro-vibrational level of the ground state is due to dipole 
coupling to all other ro-vibronic states of the ground and excited potentials. In contrast 
to the calculations in Sections [IT] and II III we do not include multi-channel effects due to 
either the hyperfine, Zeeman, or spin-orbit interaction. Instead, we base the calculation 
on our relativistic configuration-interaction MOL-RAS-CI determination of adiabatic n^fl^ 
potentials and relativistic transition dipole moments d(R) between the ground- and excited 
states. The relativistic configuration-interaction theory treats the spin-orbit interaction non- 
pertubatively for the electronic wavefunction. Consequently, we use for the polarizability 

a(hw) = j^-— E \{^v'J'M'\d{R)R-t\nvJM)\ 2 (3) 
47re c n , v ,j, M , 

f 1 11 

X I ; jz: r= + 



1^ Eni v iji — i r Yn' v 'j>/2 — (Eq vJ + hui) (.Efjw — i r Yn> v >j'/2 + tko) — Eq vJ J 

where R is the orientation of the interatomic axis, \QvJM) and \Q'v'J'M') are the ro- 
vibrational wavefunctions of initial Q and final Q' states, respectively. Here, M and M' 
are the projections of J and J' along a laboratory fixed axis. The vector e is the polar- 
ization of the laser, Eq v j is the ro-vibrational energy in the ground Q state and En> v >j' 
is the ro-vibrational energy of the excited Q' states. Contributions from scattering states 
or continuum of the excited Q' states are also included. The widths 7nvj' describe the 
spontaneous decay rate. 

In our calculation of the polarizability dipole transitions to ro-vibrational levels within 
the a 3 E + or potentials as well as to ro-vibrational levels of excited 2(0 + /~), 3(0 + /~), 

4(0+/-), 5(0+/-), 2(1), 3(1), 4(1), 5(1) and 6(1) potentials are included. The 0" states do 
not contribute to the polarizability of the X X S + ro-vibrational levels. 

If a molecule is in a ro-vibrational level of the ground electronic potential and the fre- 
quency of a laser generated optical trap is nearly resonant to some molecular transition, this 
will lead to transfer of population to the ro-vibrational level of an excited potential, which 
then by the spontaneous emission can decay to many ro-vibrational levels of the ground 
potential. As result, we lose control over the molecule in the trap. To avoid this we must 
select trap frequency intervals in which resonant excitation is unlikely. We focus on the 
dynamic polarizability of the initial and finals levels relevant to the Raman transition shown 
in the right panel of Fig. [TJ Most experimentally used trapping laser use frequencies that 
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Laser frequency (cm ) 

FIG. 8: Real (panel a) and minus the imaginary (panel b) part of the dynamic polarizability a of 
the v=— 1, —2, and —3, J = ro-vibrational levels of the a 3 S + ground state of KRb as a function 
of laser frequency. The most-widely-used trapping laser wavelengths, Ai = 1090 nm, A2 = 1064 
nm, and A3 = 1030 nm, are indicated. 

are located in the range between 9000 cm -1 and 9800 cm" 1 . In particular, we examine the 
polarizability at laser wavelengths of 1090 nm, 1064 nm, and 1030 nm. 

Figure [8] shows the real and imaginary part of the polarizability of the last three J = 
ro-vibrational levels (v = —1, —2, —3) of the a 3 S + potential as a function of laser frequency. 
The polarizability has only been evaluated every 0.05 cm" 1 . Most of the resonances in 
Fig. M are due to vibrational levels of the 3(1) potential. In the figure we have assigned the 
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FIG. 9: Real part of the dynamic polarizability a of the v=0 J = level of the X E + ground 
state of KRb as a function of laser frequency within the range of 9100 cm -1 to 9710 cm -1 . 

resonances due to transitions from the v — — 1 level of the a 3 £ + potential to even or odd 
vibrational levels of the 3(1) excited potential. The even v levels have a stronger dipole 
moment and thus wider resonance. Away from the resonances the imaginary part of the 
polarizability is six to seven orders of magnitude smaller than the real part. This indicates 
that loss due to spontaneous emission from the excited state is negligible. We also indicate 
three trap wavelengths. As we can see from Fig.[8j the wavelength Ai = 1090 nm, used in the 
JILA experiments [TJ E], is far away from the resonances and has a very-small decoherence 
rate. The two other trap wavelengths, A2 = 1064 nm and A3 = 1030 nm, lie around stronger 
resonances and their loss rate is predicted to be ten to hundred times larger. 

Figure [9] shows the real part of the polarizability of the v =0 J = level of the X 1 E + 
state of KRb as a function of laser frequency. Only one resonance is visible. It is due to 
the lowest ro-vibrational level of the excited 1 3 II(0 + ) potential. This contribution to the 
polarizability is due to the spin-orbit mixing with the 2 1 S(0 + ) potential. The benchmark 
with the shortest wavelength lies very close to this resonance. However, the precise location 
of the resonance is unknown. 

In summary, we performed a spin-coupling treatment that describes the hyperfme and 
Zeeman structure of the most weakly-bound and the most deeply-bound vibrational levels 
of the ground X 1 S + and a 3 E + potentials. The agreement with the observed experimen- 
tal structures [2] is very good. Detailed structural understanding and assignment of the 
vibrationally-cold molecules is key for the coherent control of their interactions. In ad- 
dition, we used a multi-channel description of the excited vibrational levels that includes 
independent spin-orbit coupling between multiple non-relativistic 2S+1 A ± potentials. The 
spin-orbit coupling constants were obtained from an electronic structure calculation. Finally, 
we examined the dynamic polarizability a of vibrationally cold KRb molecules as a function 
of laser frequency u. Based on this knowledge, laser frequencies can be selected to minimize 
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decoherence from loss of molecules due to spontaneous or laser-induced transitions. 
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